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Several methods have been devised to measure the weak phase 7 using decays of the type 
B — > DK ± , where it is assumed that there is no mixing in the D° — D° system. However, 
when using these methods to uncover new physics, one must entertain the real possibility that the 
measurements are affected by new physics effects in the D° — D° system. We show that even values 
of id and/or yo around 10 -2 can have a significant impact in the measurement of sin 2 7. We discuss 
the errors incurred in neglecting this effect, how the effect can be checked, and how to include it in 
. . . the analysis. 
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<D 1 I. INTRODUCTION 

Q: 

ly-j I In the next few years, the SM description of the charged current interactions through the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix Q, and, in particular, the nature of CP violation, will be subject to tests of unprecedented 

t— I ' precision. The final objective is to over-constrain the CKM matrix and, thus, probe any effects due to new physics. 
Two important tests will be the determination of sin(2/3) from the CP violating asymmetry in Bd — > J/ipKs, where 

^ ; f3 = arg {-V cd V* b /Vt d V t t), and the search for B° s - mixing. 

, Interesting constraints will also arise from the decay chains B ± — ► DK^ — > fuK ± . The idea here is based on the 

■ fact that 

OV -> D°K +) _ V: b V cs a 2 

OV A{B+^D°K+) V* b V us ai 

■ 

O ' If the final state fn is common to D° and D°, then these two amplitudes interfere and one probes essentially the 

weak phase 7 = arg (— V u dV* b /V c dV* b ). The phenomenological factor |o,2/ai| ~ 0.26 [^| accounts for the fact that the 
' B + — ► D°K + decay is color suppressed, while the B + — > D°K + decay is not. 

Inspired by a triangular construction due to Gronau and London |3j, Gronau and Wyler (GW) proposed a method 
to extract sin 2 7 which uses the decay chains i? ± — > DK^ — > fcpK^, where f cp is a CP eigenstate j|. Atwood, 
Dunietz, and Soni (ADS) have modified this method, using only final states fn which are not CP eigenstates 
Recently, Soffer has stressed the experimental advantages of combining the two strategies into a single analysis, while 
pointing out the complications due to the discrete ambiguities inherent in these methods and other measurements of 
direct CP violation [|. 

The nice feature of these decays is that they involve only tree level diagrams and, thus, are not subject to penguin 
pollution. However, one must consider what effects the mixing in the D° — D° system might have on the i? ± — > 
DK^ — > JdK^ decay chains 0, especially if these measurements are used to uncover new physics. Otherwise, new 
physics effects in the D° — D° system could be misidentified as new physics in the Bd system. In fact, Meca and Silva 
have used xd ~ 10~ 2 to argue that these effects could be of order 10% in B^ — > JdK^ decays, and they may be 
as large as 100% in B ± — * JdP ± or B ± — * decays [||. 

The main objective of our article is to study in detail the effect of D° — D° mixing on the measurement of 7 in 
S ± — > DX ± decays. If D° — D° mixing is observed, it must be incorporated into the B^ — > DX^ analysis. If only 
upper limits on mixing are known, their effect should be included as a systematic error. We also mention briefly how 
the effect of D — D° mixing may, under certain conditions, be detected in the measurement of 7. 

In section [l| we establish our notation. In section III we present the complete expressions for the B^ — > DX^ 



JdX^ decay rates in the presence of D° — D° mixing. In section IV, we start by studying the influence of D° — D° 
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mixing on the GW and ADS methods separately, concentrating on some regions of parameter space. We show that 
these effects depend on the specific value of 7 and that they might affect the extraction of sin 2 7 by as much as 
75%, even for values of xd and/or yn around 10~ 2 . Then, we combine the GW and ADS methods into a realistic 
experimental analysis, performing a scan over parameter space to discuss the impact that the mixing effects have 
on such experiments. We also show how to include the mixing effects in the analysis. We draw our conclusions in 
section |v| For completeness, the formulae relevant for the study of the D° — D° system are included in appendix JaL 
Appendix ^ contains a comparison between CP-even and CP-odd corrections to the extraction of sin 2 7. Appendix]^ 
presents an analysis of the measurements of the strong phases in the D decays that enter in the extraction of 7, and 
which may be performed at the tau-charm factories. 



II. ASSUMPTIONS AND NOTATION 



A. Parametrizations of the decay amplitudes 

One might be surprised by the fact that Eq. (|I|), is not invariant under the rephasing of the u and c quarks. In 
fact, the ratio measured experimentally in the GW and ADS methods is rather 

A(B+ -► D°K+) A(D° ->f D )' 

This ratio depends on the weak phase in Eq. (^), on the relative weak phase between the decay amplitudes A(D° — > /d) 
and A(D° — ► fr>), and it has the correct rephasing-invariant properties. The weak phase in Eq. (Q) is essentially 
given by 7. Indeed, tree level, VK-mediated D decays only probe the weak phase in the first two families, x' — 
ar g ( — V us V* d /V cs V* d ), and this is also the weak phase that appears with 7 in Eq. (Q). In the SM, x' ues around 0.003 
radians B and its presence is completely irrelevant. New physics could, in principle, affect this result by altering \' 
or by allowing for new diagrams to drive the D decays. However, both x' llOj and any additional contributions to D 
decays |If| ] are likely to remain small in most models of new physics. We will neglect them henceforth. 
For simplicity, we will use the parametrizations 

A{B+ -> WK+) = B, A{B+ -» D a K+) = eBe^e'^, 

A(B~ -» D°K-) = B, A{B- -^D°K~) = e B e' 1 " 1 e iA ° , (3) 

for the initial B+ decays; 

A(W -» K+n-) = A(D° -» K-n + ) = A, 

A{D° -> K+n-) = A(D° -» R-tt + ) = -eAe tA °, (4) 

for the D decays into non-CP eigenstates; and 

A(D° -> / cp ) - A cp , A(Tfi — > / cp ) = r)fA cp , (5) 

for D decays into a CP eigenstate f cp with CP eigenvalue r/f. In these parametrizations, we have removed all 
irrelevant overall CP-even phases. However, the differences between CP-even phases in competing paths, Ab and 
Ad, are physically meaningful, e and e are discussed in the next section. 



B. Estimates of the parameters 

In the SM, the value of 7 is already constrained by an overall analysis of the unitarity triangle. In recent reviews 
Ali and London find (l2| 

36° < 7 < 97° =>• 0.35 < sin 2 7 < 1.00 , (6) 

while Buras Jl^] quotes 44° < 7 < 97°. The variations found in the literature are mostly due to different estimates 
for the theoretical errors and to the different methods used to combine the theoretical and experimental errors. For 
definiteness, we will use Eq. ^) as a reasonable estimate for the allowed region in the SM. However, we stress that 7 
is allowed to take any value in our analysis. We are concerned not only with the effects that D° — D° mixing may 
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have on SM values of 7, but also with the possibility that such mixing effects may 'hide' new-physics by bringing 
values of 7 outside the SM allowed region into this region. 

Using iKb/Vcbl ~ 0.08 0, Eq. (Q), and assuming factorization, we get 

e ~ \V: b V cs /V c * b V us \ loa/oxl ~ 0.09 . (7) 

Notice that the parametrization in Eq. (||) is valid for any B + — ► DX + decay, with X + = K + ,tt + , p + , and higher 
resonances, but the exact values of B, e and As will vary from one channel to the next. For example, for the 
B + — > Dir + and B + — » Dp + decays we expect 

^ ~ \V; b V cd /V c * b V ud \ \a 2 / ai \ ~ 0.09 x (0.22) 2 ~ 0.004. (8) 



In addition, B[B + -> D°p+] = (1.34 ± 0.18) x 10~ 2 , B[B+ -► D a n+] = (5.3 ± 0.5) x 10~ 3 (l5|, and £[5+ 
D°is: + ]/6[S+ -> Z>°7r + ] = 0.055 ± 0.015 @. While the factor 5 2 is roughly 18 (46) times larger in B+ -> L>7r+ 
(_B + — > Dp + ) decays than it is in B + — » DK + decays, the decay chains of interest to us will scale instead with the 
B 2 e 2 factors in the CKM-suppressed decays 

B[B+ -> L>°^+] - £ 2 B[B+ -► Z?o^ + ] - 2.4 x 10~ 6 , 
-> £°7r+] - e 2 -> ^Tr + ] ~ 1.0 x 10" 7 , 

2?[£+ -> - 4 -» ^°P + ] ~ 2 - 5 x 10_7 - ( 9 ) 

As a result, the J3 + — > DK + decays are the best to extract sin 2 7. 

Similarly for the decays used in the ADS method, e is the magnitude of the ratio of the amplitudes of the doubly 
Cabibbo suppressed (DCS) decay D° — » K~tt + to the Cabibbo allowed (CA) decay D° — > K~w + . This may be 
estimated from jl7| 



^ f^ZK-^ ] ~ Vmm ~ 0m - <10) 

The parametrization in Eq. (||) is valid for any D — > /d decay, but the exact values of ^4, e and A^> will vary from one 
channel to the next. In particular, Eq. (Q) reduces to Eq. (g) when Ac = and — e = rjf = ±1 is the CP eigenvalue 
of / cp . For the CP-even eigenstates, such as K + K~ and ir + ir~ , e = — 1. For the CP-odd eigenstates, such as Kstt 

and Ks<j>, e = +1. 

The mixing in the D° — D° system may be parametrized by the mass difference divided by the average width 
(xr> = Am/r), by the width difference divided by twice the average width (yr> = Ar/(2r)), and by a CP-violating 
phase (0d) defined by qd/pd — e 2z0D . Assuming 60 — 0, 

xl+y 2 n< (6.7 x HT 2 ) 2 , (11) 

at the 95% C. L. jisj| . This bound is likely to remain stable even when one allows for Op ^ |is|j . 

Predictions for xd and yn within the SM vary considerably among the different authors pc| , but it is probably safe 
to estimate an upper bound around a few times 10 -3 . This uncertainty is due to the role played by SU(3) breaking 
effects in the long-distance part of the calculation. For example, Bucella, Lusignoli and Pugliese |H| have estimated 
the SM value for yv to lie around 1.5 x 10~ 3 . However, this result might be subject to sizeable errors, since it comes 
about through a large cancellation between two individual contributions, each of order 3 x 10 -2 . On the other hand, 
in the SM, the value of 9u is related to x' and is guaranteed to be extremely small. 

When one goes beyond the SM, one finds many models for which Op may be large and xd ~ 10 -2 , while yo 
(which is closely related to the decay rates, where one would hardly expect any large new physics contributions) 
is likely to retain its (rather uncertain) SM value PQ] . Ultimately, these values will be determined experimentally 
at _B-factories, in fixed target experiments, and at tau-charm factories. We take the point of view that, until a 
determination of xd and yp is available, their upper bounds must be included as a systematic uncertainty in the 
experimental determination of 7 from the B^ — > DK^ decays. 

The relevant point about the notation introduced here is that e and e are of order 10 _1 , while we will take 
x d ~ Vd ~ 10~ 2 as an illustrative example. (This is of the order of the sensitivity expected in the near future p2|.) 
Therefore, any effect proportional to xd/c, JAd/S, xjj/e, or yu/e is naturally of order 10% [Q, and might be larger 
depending on the exact value of the other parameters in the problem. 

We should stress that the D° — D° mixing effects might be of order unity, or even dominate, in the B ± — ► Dtt^ 1 
and B ± — > Dp ± decays, because e w ~ 0.004 ||. As a result, even SM values around xd ~ Vd ~ 10~ 3 , would lead to 
a 10% effect in these channels |w|. 
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III. THE B ± -» DK ± -» FoiT* DECAY RATES 




D°K + 



-o + 
D K 




(kV)k + 



D K 



FIG. 1. The B ± -> 
just needs to take e = ±1. 



JdK decay chain, specified for the case of fn = K 7r + . For a CP eigenstate fn = /cp one 



The -B + — > DK + — > JdK + decay chain is shown in Fig. [ft The solid lines refer to decays and the dashed lines to 
the time evolution of the D° — D° system. The functions g+(t) and g~(t) are discussed in appendix [X| and describe the 
flavor preserving and flavor changing time evolutions, respectively. The corresponding decay amplitude is obtained 
simply by adding the four possible decay paths 
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(12) 



The magnitude squared of this expression yields the time dependent decay rate. 
Using the notation of section O, we find, 



— T[B+^DK+^f D K+] = 

\g+(t)\ 2 [i 2 + e 2 - 2«cos( 7 + A B - A D )] 
+ \g-{t)\ 2 [l + e 2 ? 2 -2eecos(7 + 4^ + A s + A D )] 
+ 2 Im (<?;(%-(*)) [e(l-e 2 )sin( 7 + 2^ + A B ) - e(l ■ 
+ 2 Re (gl(t)g- (t)) [e(l + e 2 ) cos( 7 + 26 D + A B ) - e(l 



e~ 2 ) 



c 2 



sin(26»£, - 
cos(2^d 



Ad)] 
- A^)] . 



(13) 



To obtain the expression for the CP-conjugated decay rate, B — > Z?Ji" — > /u-ftT , we simply substitute 7 — > —7 
and — * — By). As shown in appendix |Aj, the time integrated decay rates may be obtained through the substitutions 
|<7±WI 2 - G± and g%{t)g-{t) - G+_. 

These expressions are completely general; no approximations were made (except for the use of \q/p\ = 1). In 
subsequent derivations we will often simplify the expressions, using the fact that xd, Ud, and e are small when 
compared with one. However the plots and estimates presented in this article were calculated using the complete 
formulae. We will also expand in e, except for decays into CP-even (CP-odd) eigenstates, where e = —1 (e = +1), 
A D = 0, and A -> A cp . 

Eq. ( |l3| ) exhibits the usual three types of CP-violating terms.F] The first line on the RHS contains a term proportional 
to sin 7 sin As, which is due to direct CP violation in the B^ — > DK ± decays (and appears multiplied by cos Ad 
in the full decay chain). The last line contains a term proportional to sin (29 d) sin Ac, which is due to direct CP 
violation in the D — > fjj decays. In both cases, CP violation requires a non-zero strong phase. The usual CP violation 
in D° — D° mixing has already been neglected in Eq. (|T^ ) due to the assumption that {qd/pdI = 1- The third line on 
the RHS contains a term proportional to sin (29jj) cos A B . This term is due to the interference between the mixing in 
the D° — D° system and its subsequent decay into the final state fjj. This term is not zero even if the strong phase 



1 These terms are better highlighted by constructing the usual CP asymmetry. But they are also present in Eq. ([L3|), along 
with CP-conserving terms. 
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vanishes, but it requires a nontrivial new phase in the mixing, 9]j. We could name this a first-mix-then-decay 
type of interference CP violation. 

In addition, Eq. ( |l3| ) contains a fourth type of CP-violating term. This appears on the third line of Eq. ([H]) as the 
term involving sin (7 + 26 d) cos Ab ■ This term persists even if As = (meaning that no strong phase is required), 
or/and 8jj = (meaning that no new phase is required in D° — D° mixing). In fact, in the limit 9d = Ad — 0, this 
term is proportional to sin 7 and, therefore, it is large even within the SM. This fourth type of CP violation was first 
pointed out by Meca and Silva , and it arises from the interference between the B — > D decays and the subsequent 
D° — D° mixing. We might call this the first-decay-then-mix type of interference CP violation. Its effect is best 
highlighted by considering a final state fjj which tags the flavor of the D° meson (corresponding to e = in Fig. [l]). 
In that case Fig. ^ has only two paths and there is still CP violation proportional to sin 7, meaning that it can be 
large even within the SM. This example should help in clearing some frequent misconceptions found in the literature. 
We see that 

• one may in principle use a single charged B decay (and its CP conjugate) to probe a source of CP violation 
which does not require strong phases, as long as a neutral D meson system is present in the intermediate state; 

• one can have a large CP-violating phase involving the D° — D° system, even within the SM. It is true that, in 
the SM, the CP-violating phases present in the D° — D° mixing and D decays are very small, because they are 
proportional to \' ■ However, there is a large CP-violating phase in the first-decay-then-mix type of interference 
CP violation involved in B — ► D decays, even within the SM; 

• the first-decay-then-mix type of interference CP violation can, in principle, be probed even when the D meson 
is detected through its (flavor-tagging) semileptonic decay. 

In general these effects require values of xd and/or yjj around 10~ 2 , in order to have an impact on i? ± — > D°K ± 
decays. Nevertheless, as we shall see below, some values of the parameters are possible for which yu ~ 10~ 2 would 
give a 75% effect, meaning that yp ~ 10~ 3 would still give a 7% effect. 



IV. THE DETERMINATION OF 7 IN THE PRESENCE OF D° - D° MIXING 

A. The no-mixing approximation 

If there were no mixing in the D° — D° system, then we would have g+(t) = e _rt//2 and <?-(t) = 0, cf. appendix [a|. 
This would leave only the uppermost and lowermost (unmixed) decay paths of Fig. [I]; the mixed paths represented by 
the diagonal dashed time-evolution lines would be absent. This is precisely what one assumes in both the GW and 
ADS methods. In that case, the decay amplitude reduces to the first two line of Eq. (jl^) and the only relevant phase 
is that in Eq. (||). As discussed above, the corresponding weak phase is essentially given by 7. 

Obviously, under the no-mixing approximation, these decays are completely insensitive to any CP-violating phase 
6d that might be present in Z?° — D° mixing. 

In the GW method, one uses time-integrated decays rates into CP-eigenstates, given by 



r [B+ -> f cp K+] <xl + i 2 + 2gcos( 7 + A B ), 
T [B- -> f cp K~] ex 1 + e 2 + 2gcos( 7 - A B ). 

Gronau and Wyler j| assume that A cp , B, and e are known from the decay rates of D 
B + — > D°K + , respectively. Therefore, Eqs. (|l4|) determine 



/ cp , B+ -> DOjfH 



(14) 
and 



c + = cos(7 + As), and 



cos(7 - A B ), 



(15) 



from which one may extract, 



sin 2 7 



sm 



2 A 



1 - c+c_ ± 



l- C+ c_ T W(l-4)(l- 



(16) 
(17) 



As discussed below, this determines 7 up to an eight-fold ambiguity 
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Unfortunately, the GW method is difficult to implement for two main reasons. The first reason is due to the 
hierarchy between the two interfering amplitudes presented in Eqs. (|l|) and (Q). This is easily seen by noting that the 
GW method hinges on extracting an interference of order e from an overall rate of order one, as shown in Eqs. (|Ti|). 
Since Eqs. ( |l4|) can be visualized as two triangles in the complex amplitude plane, this problem is sometimes explained 
by pointing out that the two triangles are squashed. 

The second reason arises from the fact that it is very difficult to measure B[B + — > D°K + ] (and, thus, e) ex- 
perimentally. One could envision determining B[B + — ► D°K + ] through the unmixed semileptonic decay chain, 
B + — > D°K + — > {X~1 + vi)dK + . However, since full reconstruction is impossible in semileptonic decays, this process 
is subject to daunting combinatoric backgrounds. Therefore, one is led to probe B[B + — > D°K + ] through the subse- 
quent decay of the D° into hadronic final states, such as K~ir + (mr) a . However, as pointed out by Atwood, Dunietz 
and Soni ||, the two (unmixed) decay chains, B+ -> D°K+ -> (K-tt+) d K+ and B+ -> WK+ -► (K-w+) D K + 
interfere at order one, since 

A(B+^D°K+)A(D°_^K-n+) ^ f ^ 1 g (lg) 
A(B+ -> D°K+)A{D° -> K-7T+) e 

This means that it is very difficult to determine e and, thus, to implement the GW method. 

Atwood, Dunietz and Soni have turned this order one interference problem into an asset. In their method, one 
uses two final states fo for which D° — > fn is a DCS decay, while D° — > Jd is Cabibbo allowed. Examples include 
fo = K~n + , if _ 7r + 7r°, i4'~7r + 7r~7r + , etc. Then, there are four decay rates, 

T [B+ « ? 2 + e? - 2«, cos( 7 + A B - A Di ), 

T [B- -» /DiJC-] « e~ 2 + e 2 - 2« 2 cos( 7 - A B + A Di ), (19) 

which may be solved for the four unknowns: e, 7, A 5 — Ab>i, and A B — A^- For each final state, i = 1 or i = 2, 
one can use the analogue of Eq. ( |l6| ) to obtain sin 2 7 and sin 2 (As — Ad%)- These expressions, of course, depend on 
the unknown e, which is determined (up to discrete ambiguities) by requiring that the expressions for sin 2 7 found for 
i = 1 and for i — 2 match. Here, although the interference terms contribute at order ee, the other terms are e 2 , and 
e 2 . As a result, the interference is of order one, and the corresponding triangles in the complex amplitude plane are 
no longer squashed. 

Soffer Q has proposed to maximize the analyzing power of the analysis by combining the GW and ADS methods, 
allowing each of them to contribute the information it is most sensitive to. In this scheme, one measures e, A^, Ab, 
and 7 by minimizing the function 

+CPConj. (20) 

X 2 compares the measured integrated decay rates in both the GW and ADS methods, with their theoretical 
expectations, T^. are functions of e, Ab, 7, and Ad, as given by Eqs. jl^ ) (for m — GW) and ( |l9| ) (for 

m = ADS). The comparison is done with respect to the measurement uncertainties, a^ p '. 

This analysis has several advantages over the individual GW and ADS methods, and, therefore, it is most likely to 
be used in the actual experiment. First, combining the relatively high-statistics GW modes with the small ADS decay 
rates improves the measurement sensitivity, due to the addition of independent data. Second, a single ADS mode is 
enough to measure all four unknowns, leaving other modes to add redundancy and statistics. Third, this analysis is 
useful in removing some discrete ambiguities, as discussed below. 



B. Discrete ambiguities 

We will now discuss the discrete ambiguities involved in the determination of 7. We start by recalling that Eq. (|l6| ) 
determines 7 up to an eight-fold ambiguity ||. A two- fold ambiguity arises from the fact that we know the signs of 
the cosines of 7 + Ab and 7 — Ab, but not the signs of the corresponding sine functions. Physically, this amounts 
to a confusion between 7 and Ab ■ A further four- fold ambiguity arises from the fact that a determination of sin 2 7 
only determines the angle 7 up to the four fold ambiguity 7, tt — 7, n + 7, and 2n — 7. Another way to interpret this 
ambiguity is to notice that Eqs. (H4) are invariant under the three independent transformations w] 





7 - 


-> -7, 


Ab — 


-A B ; 


^exchange • 


7 <- 


y A_b; 






Sir ■ 


7 - 


-> 7 + 7T, 


Ab - 


■* A B -tt 



G 



These discrete ambiguities are present in both the GW and ADS decay rates pj. The S eX change symmetry is 7 <-> Ab 
for the GW decay rates, but 7 <-> Ab — Ad for the ADS rates. Thus, by combining the two methods, the Exchange 
ambiguity can be resolved with a single ADS mode, for which Ad is large enough. Large values of Ad are expected 
in the SM. By contrast, resolving the Exchange ambiguity in the GW method, requires that A^ vary significantly 
from one B decay mode to the other. This is unlikely, because the experimental limits on CP-conserving phases in 
B — ► Dtt, D*tt, Dp and D*p Q suggest that the A B are small. 

The S'sign and S^ ambiguities are likely to degrade the value of the measurement of 7 in a non-discrete manner. 
This is due to the fact that, for 7 within the currently allowed region, Eq. @, S , s i gn S' 7r 7 tends to be close enough to 
7 to make the two values indistinguishable within the experimental errors. Overall, this results in a broad dip in the 
X 2 of Eq. ( |20|) as a function of the measured value of 7, which is likely to substantially increase the measurement 
uncertainty B . 

It is therefore instructive to notice that the term sin(7 + 29 d + Ab) of Eq. ( |l3| ) breaks the S'sign symmetry. As a 
result, one could naively expect that, incorporating mixing into the analysis would resolve the S'sign ambiguity and 
hence the SsignSV ambiguity, bringing about significant improvement over the no-mixing case. 

Unfortunately, the S s i gn -breaking term is too small to have a significant effect with foreseeable data sample sizes. 
This term vanishes in the large GW modes, in which Jd = f cp , and it accounts for about 10% of the rate in the ADS 
modes when xd = 0.01. It is estimated || that with 600 fb^ 1 collected at a next-generation B factory, the number of 
events in the ADS modes may be as large as 1300. The S'sign violating term thus accounts for some 13 events, and its 
contribution to the x 2 is only about (13/V130) 2 = 1.3, even when background is neglected. Therefore, we conclude 
that 60 times more data arc needed to effectively resolve the S S i gn ambiguity at the \ 2 — 10 level. 

C. The GW decay rates in the presence of mixing 

In the presence of D° — D° mixing, the B + — > DK + — > f cp K + decay rates involved in the GW method are altered 
according to Eq. ( |l3| ) . We find, 

r [B+ f C pK+] ex 1 + i 2 + 2r] f £cos( 7 + A B ) 

-rifX D sm(29 D ) - Vd [2ecos(7 + 29 D + A B ) + Vf cos(2# D )] , 
T [B~ -» f C pK~] <x 1 + ? 2 + 2r lf ecos( 7 - A B ) 

+T) f x D sin(2(9 D ) - y D [21 cos (7 + 29 D - A B ) + r lf cos [29 D )\ , (22) 

to linear order in xu and i/d- Recall that rjf = ±1 is the CP eigenvalue of the CP eigenstate / cp . 

In these expressions we have effectively dropped the last line of Eq. (fl2l), because the corresponding decay chain is 
suppressed both by e and by D° — D° mixing. Moreover, since we are only keeping linear terms in xd and j/c, the 
second line of Eq. (Jl^) does not contribute. These approximations were made in order to simplify the expressions. 
We stress that we have used the full expressions in our analysis and in all the computer simulations. 

The crucial step in the GW method is the identification of the interference terms on the first lines of Eqs. j2^): 
cos(7 ± Ab)- These are identified with c± in Eqs. (fl5|), which are then used in Eq. (|l6|). In the GW method 
this interference is of order e (although it is buried in an overall decay rate of order one) , while the leading mixing 
contributions are of order xd and i/d- Taking xd ~ VD ~ 10~ 2 , we expect the corrections to the determination of 
sin 2 7 to be of order xd/£ ~ Vo/e ~ 10 _1 . 

It is clear that the importance of the mixing terms depends on the exact values of 7 and A^. Therefore, they could 
be much larger than the previous naive estimate might lead us to believe. Similar considerations apply to the ADS 
method. As long as only upper bounds on D° — D° mixing are known, these effects constitute systematic uncertainties 
that must be added to any other experimental uncertainties. 

If there is D° — D° mixing, but this is ignored in the experimental analysis, then the c± extracted from the GW 
method become 

c +w = cos( 7 + A B ) - sm(29 D ) - || [2r] f ecos (7 + 29 D + A B ) + cos (29 D )\ , 

c- w = cos( 7 - A B ) + H sm(29 D ) - || \2lr) S cos (7 + 29 D -A B )+ cos {29 D )] , (23) 



2 This is the case of maximal CP-asymmetry, in which the CP-conjugate decay rate vanishes. 
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to linear order in xd and yo- The Xn terms are only important if there is also a large new CP-violating phase 9d 
in the D° — D° mixing. In contrast, the second yo term is important even if 9d is small. This is closely related 
to the fact that the linear term in the time-dependent expression for the direct D — > / cp decays is proportional to 
yo cos (29 d) ± xd sin(20D) pJJ. When Eqs. (^3|) are substituted for the c+ and c_ in Eq. (|l6|), one obtains incorrect 
values for sin 2 where ^ w stands for the 'wrong' value of 7. 

Notice that the CP-odd term proportional to xd does not involve A. However, as shown in appendix ]b|, when 
Eqs. (23) are substituted into Eq. ([l6|), the CP-odd contributions induce an error in the determination of sin 7 which 
is always proportional to sin As ■ As a result, the corrections due to yu tend to be more important than those due to 
Xd, in the small Ab limit. 

These properties are illustrated in Figs. |[ |[ and |], where we probe the effects due to xu 7^ 0, yo 7^ 0, and a 
combination of both, respectively. We have taken A# = 16.9°, which was chosen to allow a comparison with the 
results in Ref. H, and we have used the definition 



AS = sin 2 7^ — sin 2 7 



(24) 



for the difference between the 'wrong' value and the 'correct' value of sin 2 7. 
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FIG. 2. Error induced by the mixing on the GW extraction of sin 2 7. Here we take e = 0.09, As = 16.9°, 77/ = — 1, 
x D = 0.01, ijd = 0, and 6 D = 30°. 



In Fig. ^ we illustrate the effect of xd by taking rjf — —1, xd = 0.01, yn = 0, and 9d — 30°. For a true value 
of sin 2 7 ~ 0.35, we find a correction of order —0.07 (20%). Of the two possible signs in Eq. (|l6|), the plus sign, 
corresponding to the solution plotted in the figure, gives the correct value of sin 2 7, in the limit of vanishing xjj. 
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FIG. 3. Error induced by the mixing on the GW extraction of sin 2 7. Here we take e = 0.09, As = 16.9°, 77/ = +1, x r> = 0, 
y D = 0.01, and 6 D = 0°. 

In Fig. H we illustrate the effect of yn by taking 77/ = +1, Xn = 0, yr> = 0.01, and Bo = 0°. The two curves shown 
correspond to the two possible signs in Eq. ([l6]) . The solid (dashed) curve corresponds to the plus (minus) sign. For 
each value of 7, the line closest to the horizontal axis provides us with the correct value for the corrections due to 
mixing (the other solution is included in the discrete ambiguities). The correct value for sin 2 7 is best approximated 
by taking the minus sign in Eq. jl^ ) for sin 2 7 < 0.19, while it is best approximated by taking the plus sign in Eq. ( |l6| ) 
for sin 2 7 > 0.19. Therefore, the maximal value AS = ±0.14 occurs precisely for sin 2 7 = 0.19. This is a correction of 
order 75%, with either sign. This example illustrates the possibility that a value of 7 outside the SM allowed range 
might yield a value of 7™ inside the SM allowed region, due to the D° — D° mixing effects. 
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FIG. 4. Error induced by the mixing on the GW extraction of sin 2 7. Here we take e = 0.09, As = 16.9°, 77/ = — 1, 



Xd 



0.01, y D = 0.01, and G D = 30° 



In Fig. [| we take rjf = —1, ijj = 0.01, yo = 0.01, and Op = 30°. The combination of the two effects brings the 
maximum of AS up to higher values of sin 2 7, but the overall correction gets suppressed due to a partial cancellation 
between the xd and yo terms. For example, for sin 2 7 ~ 0.40, we find a correction of order 0.02 (5%). The exact 
form of this figure depends very sensitively on the precise value of 8d- As we take 0d from 0° to 90°, the figures start 
by looking like Fig ||, and end up looking like Fig 3, because the xd (vd) effect is proportional to sin26*£> (cos2#£>). 







For some values of 6d > 90° the xd and yo effects come in with the same sign and there is no partial cancellation; 
the figures become similar to Fig [|, but with larger values of AS. 
Figs. ||, ||, and |], share some important features: 

• The solution for sin 2 7 is found by using either sign in the first of Eq. ©. In fact, this is the origin of the 
7 < — > Ab discrete ambiguity in the GW method. 

• The effects ofuo on the GW method tend to be larger that those due to xd, when As is very small. As shown 
in appendix |B|, this is the result of a sinA^ suppression imposed by the inversion procedure. However, we 
should point out that this holds only under the assumption that e has been miraculously measured somehow. 
As shown by Meca and Silva @] , if one were to measure e by tagging the D meson in the final state through its 
semileptonic decay, then the xjy effect would come into the extraction of e without any sin A b suppression, and 
would be as large as the yn effects. In any case, both effects are sizeable. 

• The mixing effects may take values of 7 which are outside (inside) the SM allowed region and yield values of j w 
which are inside (outside) that region. In the first case, the mixing effects hide the new physics. In the second 
case, they give a signal for new physics when, in reality, there is none. 

• We see from Figs. |2|, [|, and that, in the presence of the mixing corrections, not all values of 7 survive the 
square root used in the inversion procedure of Eq. (|lq). Indeed, for many values of xd, yp, and 6d there is 
a range of values of 7 for which 1 — c 2 _ < or 1 — C+< 0. The worst case occurs when both 1 — c 2 _ < and 
1 - c\ < 0. In that case, the presence of the mixing will go undetected in the inversion procedure, unless one is 
independently checking whether indeed < Cj_ < 1. 

• In all cases AS = sin 2 j w — sin 2 7 — > when 7 — » ir/2 and As is small. This effect is explained in appendix 

D. The ADS decay rates in the presence of mixing 

Using Eq. ( [l3| ) we obtain the B + — > DK + — > JdK + decay rates relevant for the ADS method in the presence of 
D° — D° mixing, as 

T [B+ f Dl K+] oc e 2 + e 2 - 2e,gcos(7 + A B - A Dl ) 

-x D [e(l - e 2 ) sin( 7 + 20 D + A B ) - <h sm(26 D + A Di )] 
-y D [ecos( 7 + 20 D + A B ) ~ <h cos(26» jD + A Di )] , 

T [B~ -> foiK-] oc e 2 + e 2 - 2e l icos(- f - A B + A Dl ) 

+x D [e(l - £■) sin( 7 + 26 D - A B ) - e< sin(26» D - A Di )] 

-y D [ecos( 7 + 26 D - A B ) ~ e< cos(20 D - A Di )] , (25) 



to linear order in xd and yo- We have used here the same approximations discussed in connection with Eqs. (J2S 
and we have also expanded in 

The crucial step in the ADS methods is the identification of the interference terms cos(7 ± Ab T Ad) on the first 



lines of Eqs. (25). These terms are of order ee, while the mixing effects are of order x^e, XD^i, yo^, and yo^i- As 
a result, taking xd ~ yD ~ 10~ 2 , a naive estimate predicts the mixing effects to perturb the extraction of sin 2 7 at 
order 10%, since this is the common estimate for Xd/£, xd/ci, Vd/^-i and j/d/ci- 
If this effect is ignored in the analysis, then the c± become 

c +w = cos( 7 + A B - Ad%) + - e 2 ) sin( 7 + 20 D + A B ) - ^| sm{26 D + A Dl ) 

+|^ cos( 7 + 26 d + A B ) - || cos (26 D + A Di ), 
c_ w = cos( 7 - A B + A Dl ) - - e 2 ) sin( 7 + 26 D - A B ) + ^| sm{26 D - Am) 

cos( 7 + 26 d - A B ) - |? cos (26 D - A Di ). (26) 



3 Except for the coefficient (1 — e?), which was kept in order to allow a clear comparison between Eqs. (|25|), where (1 — ef) ~ 1, 
and Eqs. (E2), where 1 — e 2 = and the corresponding terms are absent. 
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There are two new features in Eqs. (^5|) and (^6|), which were not present in Eqs. ( p2| ) and (^). The first is the 
presence of a term proportional to id sin(7 + 20d ± As). As a result, the xn contributions no longer require the 
presence of a new CP-violating phase in the mixing Qu- The second new feature is the presence of Am- These phases 
are expected to be large in the SM. For some specific values of the parameters, the magnitude, and even the sign, 
of Arji will dramatically enhance the mixing effects. In particular, there are now effects of xu on sin 2 ^ w which are 
not proportional to sin As, but rather to sin (A^ — Ad%) or sinA^, both of which may be large. This property is 
discussed in detail in appendix^. The end result is that both xd and dd have a similar impact on the ADS method. 

Strictly speaking, in the ADS method the mixing contributions perturb the extraction of sin 2 7 both directly, as 
discussed above, and indirectly, through their effects in the determination of e. In order to obtain a simple estimate 
and to allow comparison with the GW method, we will also assume in this section that e ~ 0.09 is given (in which 
case the ADS method would require only one final state). Following Ref. ||, we use A^ = 16.9°, = ±32.3°, 
and e ~ 0.06. The effects on AS tend to mimic those shown in Figs. |[ ||, and |], except that now the xd effects arc 
present even when 60 = 0, and that all effects are typically much larger. For example, setting xd — 0.01, %)d = 0.01, 
and On — 30°, would yield A5* ~ 0.1 for the all range 0.3 < sin 2 7 < 0.8. Here, there is no cancellation between the 
xd and yD terms. 

More interestingly, one can find values of the parameters which exhibit dramatic new features. The solid (dashed) 
curves in Figs. || and ^ correspond to the plus (minus) sign in Eq. (|l6j). For each value of 7, the line closest to the 
horizontal axis provides us with the correct value for the corrections due to mixing (the other solution is included in 
the discrete ambiguities). In Fig. [|, we have taken Xu — 0.01, yp = 0, 9r> = 30°, and = —32.3°. We see that the 
largest correction has moved into higher values of sin 2 7. Introducing a nonzero yo would move this maximum even 
closer to sin 2 7=1. 
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FIG. 5. Error induced by the mixing on the ADS extraction of sin 2 7. Here we take e = 0.09, As = 16.9°, e = 0.06, 
Ad = -32.3°, x D = 0.01, y D = 0, and D = 30°. 
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FIG. 6. Error induced by the mixing on the ADS extraction of sin 2 7. Here we take e = 0.09, As = 16.9°, e = 0.06, 
A D = -32.3°, x D =0,y D = 0.01, and B D = 30°. 

In Fig. ||, we have taken Xd — 0, Vd — 0.01, 6d — 30°, and Ajj = —32.3°. Now the corrections are large and 
positive for sin 2 7 < 0.57, they are negative for 0.57 < sin 2 7 < 0.68, and they are positive again for sin 2 7 > 0.68. 
This peculiar effect is due to the fact that A 5 is minimized in different regions of sin 2 7 by each of the two possible 
signs for sin 2 7^ in Eq.jl^). The maximum deviation is ±0.1 and occurs for sin 2 7 = 0.68. This maximum occurs for 
higher values of sin 2 7 and, therefore, it is proportionally smaller (15%). We can increase this maximum by changing 
xd to 0.01. In that case, the first AS > portion of Fig. || is very suppressed, and it turns negative at sin 2 7 ~ 0.22. 
The maximum deviation would occur at sin 2 7 ~ 0.63, but would be enhanced to ±0.15 (a 24% effect). 

We have not shown ADS figures with large effects (in percentage) because they do not show new features. Neverthe- 
less, such possibilities do exist. For example, if we take Ajj = +32.3°, xn = 0.01, yn = 0.01, and Op = 60°, we obtain 
a figure which is very similar to Fig. 0, except that large deviations around 0.1 extend over a much larger domain, 



going approximately from sin 2 7 
and is ±0.13 (a 70% effect). 



0.1 up to sin 7 ~ 0.6. In this case, the largest deviation occurs at sin 7 = 0.18 



E. Combining the ADS and GW methods 



To most closely simulate the actual experiment, we will now combine the ADS and GW methods, using Eq. (J2C 
We have seen above that the GW and ADS methods are affected differently by non-trivial mixing. One may therefore 
expect that when the two methods are combined, the total effect will be smaller than in the worst-case scenario for 
either method. 

To calculate j w — 7 in this scheme, we vary Ajj, Ab, 7, Qd> xd and yjj over their allowed ranges. For each set of 
input values, we calculate of Eq. ( p0| ) using the correct expression, Eq. (|l3|), but conduct the 'wrong' analysis by 
using Eqs. ( |l4| ) and ( fl9| ) to calculate Tq^, and r^ DS , respectively. Having thus neglected mixing in our analysis, we 
proceed to minimize Eq. ( p0[ ) to obtain a measurement of j w 

The resulting distributions of A7 = 7^,-7 are shown in Fig. 0. In Fig. 0b we restrict A^ to the more likely range 
I sin(As)| < 0.5. Since the actual values of the phases are not known, it is not meaningful to analyze the distributions 
of Fig. in detail. These distributions clearly indicate, however, that neglecting to account for D° — D° mixing may 
significantly bias the measurement of 7. 



4 As in j6|, the input value of e is taken from Eq. 
other observables) when we minimize Eq. (pp|). 



(Q), but its 'experimental' output value is then found (together with the 
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FIG. 7. Distributions of A7 = j w — 7, the error in the measurement of 7 due to failure to incorporate mixing into the 
analysis. Parameters are varied in the ranges —180° < Ad < 180°, 36° < 7 < 97°, 9d = 0°, 30°, or 45° xd — or 0.01, and 
y D = or 0.01. A B is varied within -180° < A D < 180° (a) or | sin(A s )| < 0.5 (b). 



As we have seen before, the measurement of 7 depends on the the precise values of the parameters in the D 
system. Given the current bounds on xd and y n , we expect the measurements of e and Ac to be rather insensitive to 
D° — D° mixing. The measured value of e, c.f. Eq. (|l0|), is already used in the ADS method. Once Ad is measured, 
c.f. appendix 0, it can also be used as a known input parameter in the fit for 7. We have stressed the fact that these 
decays also depend on id, yD, and 9d- The more we know about these quantities, the better will be our bound on 7 
from the _B ± — > DK^ decays. 



V. CONCLUSIONS 



The decay chains — > DK ± — > JdK^ provide a good opportunity to determine the CKM phase 7. Naturally, 
information about the D decays must be included in the analysis, either as parameters to be determined from the 
overall fit or as fixed quantities known from other D system experiments. This is well appreciated for the DCS decay 
parameter e and for the strong phase Ad. 

In this article, we stress the fact that this is also true for the parameters involved in D° — D° mixing. This point 
cannot be (as it often is) overlooked, because the extraction of 7 hinges on measurements of small quantities. In the 
GW method there is a small interference; in the ADS method the decay rates are small. As a result, mild values of 
xd and/or yD can have an important effect in these methods. 

We have shown that dramatic effects are indeed present when the GW and ADS methods are used individually. In 
general, combining these methods reduces the errors involved, but one does still find deviations of order j w — 7 ~ 10°. 
These effects may simulate the presence of new physics by taking values of 7 inside the SM allowed region into values 
of j w outside that region. They may also obscure the presence of new physics by taking new physics values of 7 
outside the SM allowed region and yielding j w inside that region. The importance of this error in the determination 
of 7 is made more acute by the discrete ambiguities associated with the GW and ADS methods. 

As a result 
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• the determination of 7 with the decays B — > DK must be made in connection with the search for mixing in the 
D — D° mixing; 

• any uncertainties due to lack of knowledge of xjj, yu, and must be correctly included as systematic uncertainties 
to the extraction of sin 2 7. 

• once Xd, yo, and 6 are determined experimentally, they can be included in the analysis as known parameters. 
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APPENDIX A: TIME DEPENDENT EVOLUTION OF THE D° - D° SYSTEM. 

Assuming CPT invariance, the mass eigenstates of the D° — D° system are related to the flavor eigenstates by 

\D H ) =p D \D°)+q D \W) , 

\D L )=p D \D°)-q D \D°) , (Al) 

with \pd\ 2 + \qd\ 2 = 1 and 



q D _ Am - |Ar [R^ 



■ ( A2 ) 

Pd 2i?i2 V R12 

where Am = m# — m^ (7J-heavy, L-light) is positive by definition, Ar = r# — T^, and i? 12 is the off-diagonal matrix 

element in the effective time evolution in the D° — D° space. 

Consider a D° (D°) meson which is created at time t — and denote by D°(t) (D°(t)) the state that it evolves 
into after a time t, measured in the rest frame of the meson D. One finds 

\D°(t))= g+ (t)\D°) + ^(t)\W), 

PD 

\W(t)) = PRg_(t)\D°)+g + (t)\D»), (A3) 

where 

ff± (t) = i( e -^*± e -^*), (A4) 

fJ-H = rn H - iT H /2, and fi L = m L ~ iT L /2. 

It is useful to trade the mass and width differences for the dimensionless parameters x = Am/T and y = AT/ (2T), 
where T = (Th + Fl)/2 is the average width. We already know from studies of the direct D decays, D — > K + tt~ , 
that x 2 D + y 2 D < (6.7 x 10~ 2 ) , at the 95% C. L. jl8). Therefore, we may expand the time-dependent functions as 

g+{t) ~ e- mt e- T ' 2 [l + (x D - i Vd ) 2 t 2 /A +•••], 

, 9 _ (t) ~ e~ mt e- T ' 2 l(-ix D - y D )r/2 + •••], (A5) 

where r = Tt is the (proper) time of the D evolution, in units of the D average width. 
The time-dependent decay rates involve 
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\9+(t)f 

\g-(t)\ 2 



-n 



-n 



, ^ /A X 

cosh — h cos (Ami) 



, /A X 

cosh — cos (Ami) 



[l + (y 2 D -xl)r 2 /4+---] 



g* + (t)g_(t) = -[e- VH " - < 



The time integrated decay rates involve 



t „-v L t 



2ie~ rt sin(Ami)] ~ e _r [(-£/£> - ix D )r/2 + 



G+ = J o \g+(t)\ z dt = — 



i / i 



i 



2T \ 1 - y 2 ' 1 + x 2 ) " T V 2 



/■+oo 1 / 1 

G- = ] o IM*)l a * = ^( T 



1 



1 



G + _ 



+ 00 



g* + (t)g_(t)dt 



-y 2 1 + x 2 J T V 2 

J_ f ~ y ~ ix \ ~ Vd ~ iXD 

2T \ l-y 2 l + x 2 ) ~ T 



(A6) 
(A7) 

(A8) 
(A9) 

(A10) 
(All) 
(A12) 



APPENDIX B: CP-EVEN AND CP-ODD CORRECTIONS TO THE EXTRACTION OF sin 2 7 



We have seen in Eqs. (|23|) and (|26|) that the presence of mixing may affect cos(7 ± A) by 

c +w = cos(7 + A) + 5 e + S a , 
cos(7 - A) + S e - S Q , 



c_. 



(Bl) 



where S e and S stand for CP-even and CP-odd corrections, respectively. The strong phase is given by A = Ab in 
the GW method, and by A = Ab — Adi in the ADS method. 

For example, the corrections to the GW method described in Eqs. (|23|), are 



Se 
So 



~ [2lf] f cos (7 + 20 D ) cos (A s ) + cos (28 D )} , 

XD 



2c 



- sm(26 D ) + y D rjf sin (7 + 20 D ) sin (A B ). 



(B2) 



Notice that the CP-odd term proportional to xb does not involve A. In fact, that term is due to the interference CP 
violation present in the decay D — > f cp , and it is complet ely independent of the production mechanism of the neutral 
D meson. Nevertheless, as we will now show, when Eqs. (|Bl| ) are substituted into Eq. (|l6|), the CP-odd contributions 
to the difference between the wrong value of sin 2 j w (obtained using c± w ) and the correct value of sin 2 7 (obtained 
using cos(7 ± A)), are always proportional to sin A. 

We start by noting that c± w have a standard CP-even (CP-odd) component given by cos 7 cos A (— sin 7 sin A), in 
addition to the mixing component 5 e (5 ). We will denote their sum by 



and write 



Then 



c e = cos 7 cos A + S e , 
c = — sin 7 sin A + S , 



(B3) 



(B4) 



sin 7„ 



1 - C +W C- W ± - 4J(1 - Ct w ) 

1 - c 2 + c 2 ± V(l-c 2 + c 2 ) 2 -4 C 2 



(B5) 
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and the results depend only on c 2 and c 2 .0 

Consequently, the mixing contributions to c 2 and c 2 are either quadratic (and, thus, much suppressed, although they 
could dominate for small values of cos(7 ± A)), or linear, but appearing on ly i n the combinations S e cos 7 cos A and 



(5 sin7sinA. Now, in the GW method xd only shows up in 8 , c. f. Eqs. (B2). Therefore, the mixing contribution 
to the GW extraction of sin 2 7^ proportional to ijj sin(20c) is also proportional to sin A^ and, thus, it vanishes in 
the A B -> limit. 

We can also use this discussion to explain why sin 2 7^ is very similar to sin 2 7 for 7 ~ tt/2, provided that A is 
small. Indeed, in that case c 2 = <5 2 , c 2 = sin 2 7 sin 2 A — 25 sin 7 sin A + <5 2 , and both contributions become small. 

The situation is altered in the ADS method because there we have a new CP-even phase, Am, which is expected 
to be large. Using Eqs. ( p6|) we find the leading CP-even and CP-odd contributions to c± w to be 

S e = -7^(1 - <h) co s (7 + 20b) sin A B + || cos (26 D ) sin A Dl 

S = - e?) sin (7 + 26 D ) cos A B - ^ sin (26>rj) cos A Di 

- 1^ sin (7 + 26 D ) sin A B + || cos (20 D ) cos A Dl , (B6) 

respectively. These are added to the standard contributions cos 7 cos (A^ — Aoi) and — sin 7 sin (As — Aoi), respec- 
tively. The presence of a potentially large A Di has two consequences. Firstly, the CP-odd quantity that exists even in 
the absence of mixing, is proportional to sin (A B — Aoi) and can be large. Therefore, the fact that the mixing CP-odd 
contribution to sin 2 7^ linear in 5 always appears multiplied by sin (A^ — A B i) ceases to constitute a suppression 
factor. Secondly, there is now one Xd contribution to S e which is proportional to sin A^. This will interfere with the 
standard CP-even contribution, cos 7 cos (A B — A^i), as is also unsuppressed in the A# — > limit|] 



APPENDIX C: MEASURING Ad AT A TAU-CHARM FACTORY 



3591297-001 



suppressed amplitude, Ancs, an d their interference, 




FIG. 8. Obtaining the phase A_d of Equation (CI) from the Cabibbo allowed D° decay amplitude, Ac a, the doubly Cabibbo 

Ac a ± A D cs- 



We proceed to study the measurement of Ad at a charm factory, operating at the "0(3770) resonance. For simplicity, 
we will discuss this measurement in the context of the SM. Given the current bounds 011 x B and yo, we expect this 
measurement to be rather insensitive to D — D° mixing. The relation 



Ad = arg A(D° f D )A(D° f D )* 



(CI) 



is graphically represented in Figure ||, demonstrating how to obtain A^> from the the Cabibbo allowed D decay 
amplitude, Ac a, the doubly Cabibbo suppressed amplitude, Adcs, and their interference, A± = Ac a ± Adcs- 



5 Notice that this property is completely general and holds for any method in which one is ultimately measuring the square 
of a CP-violating quantity. Indeed, the quantity sin 2 7™ is CP-even (its signs remains the same under a CP transformation) . 
Therefore, CP-even (c e ) and CP-odd (c ) contributions to this quantity cannot interfere with one another and, moreover, c 
can only contribute in the combination c 2 ,. 

6 There is also a new CP-even contribution, — e 2 ) cos (7 + 20d), but it is proportional to sin A_g. 
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While Ac a and Adcs have been measured at CLEO |L7[ for the K tt + mode by using D* + decays to tag the 
D° flavor, measuring A± requires producing D a — D° pairs in a known coherent state. It is therefore best to perform 
all three measurements at the charm factory, canceling many systematic errors in the construction of the triangles of 
Figure ||. To measure the amplitude A + (A-), one of the ?/>(3770) daughters is tagged as a D°_ by observing 

it decay into a CP-odd (CP-even) state, such as K$ it (K + K~). The other daughter is then D° + (D°_), and the 
fraction of the time that it is seen decaying into K~ir + gives the interference amplitude 4|^4±| 2 - We immediately find 



cos A_d = ± 



\A CA \ 2 - \A DCS \ 



2\A CA \\A D cs\ 



(C2) 



Due to the low statistics tagging scheme of the \A±\ measurement and the fact that I^dcs) *C |^4cm|, the error in 
cos Ajj is dominated by the \A±\ measurement error. Hence 



Ccos A £ 



a \Ai 



\A DCS \ 

1 



a \A ± \ 

\A±\ 

A CA 



A 



CA 



Adcs 



A 



DCS 



(C3) 



where Na± is the number of events observed in the A± channels, and we made use of I^cmI ~ l^±l- Since the event is 
fully reconstructed, background is expected to be small, and its contribution to u CO sA D is neglected in this discussion. 
Na ± is given by 



N A , 



N 



DD 



B{D° -> K-TT+) e(if-7r H 



(C4) 



where N D -p is the number of ■i/ , (3770) — > D°D events, ti are the CP-eigenstates used for tagging, and e(X) is the 



reconstruction efficiency of the state X. Taking || J2i 
1/V0.0031, and B(D° -> R--K+) w 0.04, Eq. (03) becomes 



U) e{ti) re 0.02, e(K-TT+) w 0.8 and \A CA /A 



DCS I 



370 



Ccos A n 



(C5) 



DD 



It is expected that in one year the charm factory will collect 10 fb _1 , or N D -p = 2.9 x 10 7 [|25|, resulting in 
a cos A D ~ 0.065. Thus, cosA^ can be measured to high precision, even in the presence of background and with 
relatively modest luminosity. We note that the same measurement technique can be used with multi-body D° decays, 
in which cos A^> varies over the available phase space. While the relative statistical error in every small region of 
phase space will be large, its effect on the measurement of 7 in B — > DK will be proportionally small. The total error 
in 7 due to AcosA^ will be as small as in the two body K~ir + mode, up to differences in D° branching fractions, 
reconstruction efficiencies, and backgrounds. 
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